Simultaneous species detection and discovery with environmental DNA metabarcoding: A freshwater mollusk case study.

Environmental DNA (eDNA) sampling is a powerful tool for rapidly characterizing biodiversity patterns for specious, cryptic taxa with incomplete taxonomies. One such group that are also of high conservation concern are North American freshwater gastropods. In particular, springsnails of the genus Pyrgulopsis (Family: Hydrobiidae) are prevalent throughout the western United States where >140 species have been described. Many of the described species are narrow endemics known from a single spring or locality, and it is believed that there are likely many additional species which have yet to be described. The distribution of these species across the landscape is of interest because habitat loss and degradation, climate change, groundwater mining, and pollution have resulted in springsnail imperilment rates as high as 92%. Determining distributions with conventional sampling methods is limited by the fact that these snails are often <5 mm in length with few distinguishing morphological characters, making them both difficult to detect and to identify. We developed an eDNA metabarcoding protocol that is both inexpensive and capable of rapid, accurate detection of all known Pyrgulopsis species. When compared with conventional collection techniques, our pipeline consistently resulted in detection at sites previously known to contain Pyrgulopsis springsnails and at a cost per site that is likely to be substantially less than the conventional sampling and individual barcoding that has been done historically. Additionally, because our method uses eDNA extracted from filtered water, it is non-destructive and suitable for the detection of endangered species where "no take" restrictions may be in effect. This effort represents both a tool which is immediately applicable to taxa of high conservation concern across western North America and a case study in the broader application of eDNA sampling for landscape assessments of cryptic taxa of conservation concern.

Patterns of Wolf Dispersal Respond to Harvest Density across an Island Complex.

Wolves are highly mobile predators and can disperse across a variety of habitats and over long distances. However, less is known about dispersal capabilities across water and among islands. The biogeography of island systems fosters spatially structured local populations, and their degree of connectivity may influence the dynamics and long-term viability of the regional population. We sought to quantify wolf dispersal rate, distance, and dispersal sex bias throughout Prince of Wales Island, a 6670 km2 island in southeast Alaska, and the surrounding islands that constitute the wildlife management unit (9025 km2). We also investigated patterns of dispersal in relation to hunting and trapping intensity and wolf population density. We used DNA data collected during 2012-2021 long-term monitoring efforts and genotyped 811 wolves, 144 of which (18%) were dispersers. Annual dispersal rates were 9-23% and had a weakly positive relationship with wolf density. Wolves dispersed 41.9 km on average (SD = 23.7 km), and males and females did not disperse at different rates. Of the dispersing wolves, 107 died, and the majority (n = 81) died before they were able to settle. The leading manner of death was trapping (97% of mortalities), and wolves tended to disperse from areas with low harvest density to areas where harvest density was relatively higher. Dispersal occurred both to and from small islands and the larger Prince of Wales Island, indicating bidirectional as opposed to asymmetrical movement, and the genetic overlap of wolf groups demonstrates connectivity throughout this naturally patchy system. Island ecosystems have different predator-prey dynamics and recolonization processes than large, intact systems due to their isolation and restricted sizes; thus, a better understanding of the degree of population connectivity including dispersal patterns among islands in the Prince of Wales archipelago could help inform the management and research strategies of these wolves.

Long-term video and genetic data yield insights into complex sociality of a solitary large carnivore.

American black bears (Ursus americanus) may be more social than currently understood. We used long-term video and genetic data to evaluate social interactions among wild, independent-aged black bear on a conservation property in western Montana, USA. We used multinomial logistic regression to evaluate predictions about male-male interactions within the context of individual fitness, female-female interactions within the context of inclusive fitness, and male-female interactions within the context of female counterstrategies to infanticide. Overall, our findings challenged the assumption that independent-aged bears interact only during the mating season or when concentrated feeding sites are present. We documented 169 interaction events by at least 66 bear pairs, 92 (54%) of which occurred outside of the peak mating season and in the absence of concentrated feeding sites. The probability that male-male pairs engaged in play and other non-agonistic behaviours was higher than that for female-female pairs. Conversely, the probability that female-female pairs engaged in chase behaviour was higher than that for male-male and male-female pairs. We documented evidence of female mate choice, female resource defense, sexually selected infanticide (SSI), and female counterstrategies to avoid SSI. Our findings improve our understanding of ursid ethology and underscore the complexity of ursid sociality.

Variation in adult and pup wolf diets at natal den sites is influenced by forest composition and configuration.

Although wolves are wide-ranging generalist carnivores throughout their life cycle, during the pup-rearing season wolf activity is focused on natal den sites where pup survival depends upon pack members provisioning food. Because prey availability is influenced by habitat quality within the home range, we investigated the relative importance of prey species for adults and pups and further examined the relationship between habitat characteristics, wolf diet, and litter size on Prince of Wales Island (POW) in Southeast Alaska. During 2012-2020, we detected 13 active den sites within the home ranges of nine wolf packs. We estimated minimum pup counts using motion-detecting cameras and individual genotypes from noninvasive samples (hair: n = 322; scat: n = 227) and quantified wolf diet composition using fecal DNA metabarcoding (n = 538). We assessed habitat composition, configuration, and connectivity within denning and annual home ranges estimated using wolf GPS-collar data. Contrary to expectations, wolves had a more constricted diet during denning season (April 15-July 31), and within this season pups had a narrower dietary niche (species richness [S] = 4) focused more on deer (relative frequency of occurrence [O/I] = 0.924) than adults (S = 15; deer O/I = 0.591). Litter size had a positive relationship with the relative frequency of deer in a wolf pack's diet. Wolf consumption of deer was positively associated with the proportion of young-growth forest (≤25 years old) within denning and annual home ranges. High levels of vegetation patch interspersion, and the density of closed logging roads were also important predictors, suggesting these habitat qualities were influential for increasing the availability of deer to wolves. Our results contrast with previous research indicating wolf pup diets included more alternate prey (i.e., beaver) than adults and emphasize the importance of deer to wolf viability on POW, especially during denning season.

Wolf Dispersal Patterns in the Italian Alps and Implications for Wildlife Diseases Spreading.

Wildlife dispersal directly influences population expansion patterns, and may have indirect effects on the spread of wildlife diseases. Despite its importance to conservation, little is known about dispersal for several species. Dispersal processes in expanding wolf (Canis lupus) populations in Europe is not well documented. Documenting the natural dispersal pattern of the expanding wolf population in the Alps might help understanding the overall population dynamics and identifying diseases that might be connected with the process. We documented 55 natural dispersal events of the expanding Italian wolf alpine population over a 20-year period through the use of non-invasive genetic sampling. We examined a 16-locus microsatellite DNA dataset of 2857 wolf samples mainly collected in the Western Alps. From this, we identified 915 individuals, recaptured 387 (42.3%) of individuals, documenting 55 dispersal events. On average, the minimum straight dispersal distance was 65.8 km (±67.7 km), from 7.7 km to 517.2 km. We discussed the potential implications for maintaining genetic diversity of the population and for wildlife diseases spreading.

Age-specific, population-level pedigree of wild black bears provides insights into reproduction, paternity, and maternal effects on offspring apparent survival.

Wildlife pedigrees provide insights into ecological and evolutionary processes. DNA obtained from noninvasively collected hair is often used to determine individual identities for pedigrees and other genetic analyses. However, detection rates associated with some noninvasive DNA studies can be relatively low, and genetic data do not provide information on individual birth year. Supplementing hair DNA stations with video cameras should increase the individual detection rate, assuming accurate identification of individuals via video data. Video data can also provide birth year information for individuals captured as young of the year, which can enrich population-level pedigrees. We placed video cameras at hair stations and combined genetic and video data to reconstruct an age-specific, population-level pedigree of wild black bears during 2010-2020. Combining individual birth year with mother-offspring relatedness, we also estimated litter size, interlitter interval, primiparity, and fecundity. We used the Cormack-Jolly-Seber model in Program Mark to evaluate the effect of maternal identity on offspring apparent survival. We compared model rankings of apparent survival and parameter estimates based on combined genetic and video data with those based on only genetic data. We observed 42 mother-offspring relationships. Of these, 21 (50%) would not have been detected had we used hair DNA alone. Moreover, video data allowed for the cub and yearling age classes to be determined. Mean annual fecundity was 0.42 (95% CI: 0.27, 0.56). Maternal identity influenced offspring apparent survival, where offspring of one mother experienced significantly lower apparent survival (0.39; SE = 0.15) than that of offspring of four other mothers (0.89-1.00; SE = 0.00-0.06). We video-documented cub abandonment by the mother whose offspring experienced low apparent survival, indicating individual behaviors (e.g., maternal care) may scale up to affect population-level parameters (e.g., cub survival). Our findings provide insights into evolutionary processes and are broadly relevant to wildlife ecology and conservation.

Molecular species delimitation refines the taxonomy of native and nonnative physinine snails in North America.

Being able to associate an organism with a scientific name is fundamental to our understanding of its conservation status, ecology, and evolutionary history. Gastropods in the subfamily Physinae have been especially troublesome to identify because morphological variation can be unrelated to interspecific differences and there have been widespread introductions of an unknown number of species, which has led to a speculative taxonomy. To resolve uncertainty about species diversity in North America, we targeted an array of single-locus species delimitation methods at publically available specimens and new specimens collected from the Snake River basin, USA to generate species hypotheses, corroborated using nuclear analyses of the newly collected specimens. A total-evidence approach delineated 18 candidate species, revealing cryptic diversity within recognized taxa and a lack of support for other named taxa. Hypotheses regarding certain local endemics were confirmed, as were widespread introductions, including of an undescribed taxon likely belonging to a separate genus in southeastern Idaho for which the closest relatives are in southeast Asia. Overall, single-locus species delimitation was an effective first step toward understanding the diversity and distribution of species in Physinae and to guiding future investigation sampling and analyses of species hypotheses.

Response of a remnant marmot population to habitat enhancement yields insights into marmot ecology.

We evaluated the response of a remnant population of yellow-bellied marmots (Marmota flaviventris) to targeted habitat enhancement in an ecological system that had been degraded during ~100 years of intensive livestock management, including marmot eradication. We used capture-recapture data and a novel use of a multistate framework to evaluate geographic expansion of the marmot population pre- and post-habitat enhancement. We also estimated age-structured survival, reproduction, and sex ratios. The marmot population appeared to respond positively to new habitat opportunities created by habitat enhancement: the number of marmots captured increased from three marmots pre-habitat enhancement to 54 (28 adults and yearlings, 26 young) post-habitat enhancement at the end of the study. Marmots expanded geographically by transitioning into habitat-enhanced areas, and adult females occupied and reproduced in all habitat-enhanced areas. The sex ratio of the young population in 2019 was strongly female-biased, which may have been influenced by poor body condition of breeding females owing to unusually prolonged snow cover that year. Adult and yearling survival were within the range of that reported for colonial adults and yearlings in Colorado. Our results suggest that active habitat enhancement can assist in the recovery of marmot populations in systems where marmots historically existed.

Pliocene-Early Pleistocene Geological Events Structure Pacific Martens (Martes caurina).

The complex topography, climate, and geological history of Western North America have shaped contemporary patterns of biodiversity and species distributions in the region. Pacific martens (Martes caurina) are distributed along the northern Pacific Coast of North America with disjunct populations found throughout the Northwestern Forested Mountains and Marine West Coast Forest ecoregions of the West Coast. Martes in this region have been classified into subspecies; however, the subspecific designation has been extensively debated. In this study, we use genomic data to delineate conservation units of Pacific marten in the Sierra-Cascade-Coastal montane belt in the western United States. We analyzed the mitochondrial genome for 94 individuals to evaluate the spatial distribution and divergence times of major lineages. We further genotyped 401 individuals at 13 microsatellite loci to investigate major patterns of population structure. Both nuclear and mitochondrial DNA suggest substantial genetic substructure concordant with historical subspecies designations. Our results revealed that the region contains 2 distinct mitochondrial lineages: a Cascades/Sierra lineage that diverged from the Cascades/coastal lineage 2.23 (1.48-3.14 mya), consistent with orogeny of the Cascade Mountain chain. Interestingly, Pacific Martes share phylogeographic patterns similar with other sympatric taxa, suggesting that the complex geological history has shaped the biota of this region. The information is critical for conservation and management efforts, and further investigation of adaptive diversity is warranted following appropriate revision of conservation management designations.

Integrative taxonomy refutes a species hypothesis: The asymmetric hybrid origin of Arsapnia arapahoe (Plecoptera, Capniidae).

Molecular tools are commonly directed at refining taxonomies and the species that constitute their fundamental units. This has been especially insightful for groups for which species hypotheses are ambiguous and have largely been based on morphological differences between certain life stages or sexes, and has added importance when taxa are a focus of conservation efforts. Here, we examine the taxonomic status of Arsapnia arapahoe, a winter stonefly in the family Capniidae that is a species of conservation concern because of its limited abundance and restricted range in northern Colorado, USA. Phylogenetic analyses of sequences of mitochondrial and nuclear genes of this and other capniid stoneflies from this region and elsewhere in western North America indicated extensive haplotype sharing, limited genetic differences, and a lack of reciprocal monophyly between A. arapahoe and the sympatric A. decepta, despite distinctive and consistent morphological differences in the sexual apparatus of males of both species. Analyses of autosomal and sex-linked single nucleotide polymorphisms detected using genotyping by sequencing indicated that all individuals of A. arapahoe consisted of F1 hybrids between female A. decepta and males of another sympatric stonefly, Capnia gracilaria. Rather than constitute a self-sustaining evolutionary lineage, A. arapahoe appears to represent the product of nonintrogressive hybridization in the limited area of syntopy between two widely distributed taxa. This offers a cautionary tale for taxonomists and conservation biologists working on the less-studied components of the global fauna.

Climate, Demography, and Zoogeography Predict Introgression Thresholds in Salmonid Hybrid Zones in Rocky Mountain Streams.

Among the many threats posed by invasions of nonnative species is introgressive hybridization, which can lead to the genomic extinction of native taxa. This phenomenon is regarded as common and perhaps inevitable among native cutthroat trout and introduced rainbow trout in western North America, despite that these taxa naturally co-occur in some locations. We conducted a synthetic analysis of 13,315 genotyped fish from 558 sites by building logistic regression models using data from geospatial stream databases and from 12 published studies of hybridization to assess whether environmental covariates could explain levels of introgression between westslope cutthroat trout and rainbow trout in the U.S. northern Rocky Mountains. A consensus model performed well (AUC, 0.78-0.86; classification success, 72-82%; 10-fold cross validation, 70-82%) and predicted that rainbow trout introgression was significantly associated with warmer water temperatures, larger streams, proximity to warmer habitats and to recent sources of rainbow trout propagules, presence within the historical range of rainbow trout, and locations further east. Assuming that water temperatures will continue to rise in response to climate change and that levels of introgression outside the historical range of rainbow trout will equilibrate with those inside that range, we applied six scenarios across a 55,234-km stream network that forecast 9.5-74.7% declines in the amount of habitat occupied by westslope cutthroat trout populations of conservation value, but not the wholesale loss of such populations. We conclude that introgression between these taxa is predictably related to environmental conditions, many of which can be manipulated to foster largely genetically intact populations of westslope cutthroat trout and help managers prioritize conservation activities.

Patterns of hybridization among cutthroat trout and rainbow trout in northern Rocky Mountain streams.

Introgressive hybridization between native and introduced species is a growing conservation concern. For native cutthroat trout and introduced rainbow trout in western North America, this process is thought to lead to the formation of hybrid swarms and the loss of monophyletic evolutionary lineages. Previous studies of this phenomenon, however, indicated that hybrid swarms were rare except when native and introduced forms of cutthroat trout co-occurred. We used a panel of 86 diagnostic, single nucleotide polymorphisms to evaluate the genetic composition of 3865 fish captured in 188 locations on 129 streams distributed across western Montana and northern Idaho. Although introgression was common and only 37% of the sites were occupied solely by parental westslope cutthroat trout, levels of hybridization were generally low. Of the 188 sites sampled, 73% contained ≤5% rainbow trout alleles and 58% had ≤1% rainbow trout alleles. Overall, 72% of specimens were nonadmixed westslope cutthroat trout, and an additional 3.5% were nonadmixed rainbow trout. Samples from seven sites met our criteria for hybrid swarms, that is, an absence of nonadmixed individuals and a random distribution of alleles within the sample; most (6/7) were associated with introgression by Yellowstone cutthroat trout. In streams with multiple sites, upstream locations exhibited less introgression than downstream locations. We conclude that although the widespread introduction of nonnative trout within the historical range of westslope cutthroat trout has increased the incidence of introgression, sites containing nonadmixed populations of this taxon are common and broadly distributed.

Cottus schitsuumsh, a new species of sculpin (Scorpaeniformes: Cottidae) in the Columbia River basin, Idaho-Montana, USA.

Fishes of the genus Cottus have long been taxonomically challenging because of morphological similarities among species and their tendency to hybridize, and a number of undescribed species may remain in this genus. We used a combination of genetic and morphological methods to delineate and describe Cottus schitsuumsh, Cedar Sculpin, a new species, from the upper Columbia River basin, Idaho-Montana, USA. Although historically confused with the Shorthead Sculpin (C. confusus), the genetic distance between C. schitsuumsh and C. confusus (4.84-6.29%) suggests these species are distant relatives. Moreover, the two species can be differentiated on the basis of lateral-line pores on the caudal peduncle, head width, and interpelvic width. Cottus schitsuumsh is also distinct from all other Cottus in this region in having a single small, skin-covered, preopercular spine. Haplotypes of mtDNA cytochrome oxidase c subunit 1 of C. schitsuumsh differed from all other members of the genus at three positions, had interspecific genetic distances typical for congeneric fishes (1.61-2.74% to nearest neighbors), and were monophyletic in maximum-likelihood trees. Microsatellite analyses confirmed these taxonomic groupings for species potentially sympatric with C. schitsuumsh and that fish used in morphological comparisons were unlikely to be introgressed. Its irregular distribution, in the Spokane River basin in Idaho and portions of the Clark Fork River basin in Montana, may have resulted from human-assisted translocation.

DNA barcoding at riverscape scales: assessing biodiversity among fishes of the genus Cottus (Teleostei) in northern Rocky Mountain streams.

There is growing interest in broad-scale biodiversity assessments that can serve as benchmarks for identifying ecological change. Genetic tools have been used for such assessments for decades, but spatial sampling considerations have largely been ignored. Here, we demonstrate how intensive sampling efforts across a large geographical scale can influence identification of taxonomic units. We used sequences of mtDNA cytochrome c oxidase subunit 1 and cytochrome b, analysed with maximum parsimony networks, maximum-likelihood trees and genetic distance thresholds, as indicators of biodiversity and species identity among the taxonomically challenging fishes of the genus Cottus in the northern Rocky Mountains, USA. Analyses of concatenated sequences from fish collected in all major watersheds of this area revealed eight groups with species-level differences that were also geographically circumscribed. Only two of these groups, however, were assigned to recognized species, and these two assignments resulted in intraspecific genetic variation (>2.0%) regarded as atypical for individual species. An incomplete inventory of individuals from throughout the geographical ranges of many species represented in public databases, as well as sample misidentification and a poorly developed taxonomy, may have hampered species assignment and discovery. We suspect that genetic assessments based on spatially robust sampling designs will reveal previously unrecognized biodiversity in many other taxa.

Historical and contemporary DNA indicate fisher decline and isolation occurred prior to the European settlement of California.

Establishing if species contractions were the result of natural phenomena or human induced landscape changes is essential for managing natural populations. Fishers (Martes pennanti) in California occur in two geographically and genetically isolated populations in the northwestern mountains and southern Sierra Nevada. Their isolation is hypothesized to have resulted from a decline in abundance and distribution associated with European settlement in the 1800s. However, there is little evidence to establish that fisher occupied the area between the two extant populations at that time. We analyzed 10 microsatellite loci from 275 contemporary and 21 historical fisher samples (1880-1920) to evaluate the demographic history of fisher in California. We did not find any evidence of a recent (post-European) bottleneck in the northwestern population. In the southern Sierra Nevada, genetic subdivision within the population strongly influenced bottleneck tests. After accounting for genetic subdivision, we found a bottleneck signal only in the northern and central portions of the southern Sierra Nevada, indicating that the southernmost tip of these mountains may have acted as a refugium for fisher during the anthropogenic changes of the late 19(th) and early 20(th) centuries. Using a coalescent-based Bayesian analysis, we detected a 90% decline in effective population size and dated the time of decline to over a thousand years ago. We hypothesize that fisher distribution in California contracted to the two current population areas pre-European settlement, and that portions of the southern Sierra Nevada subsequently experienced another more recent bottleneck post-European settlement.